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Abstract — [] In this paper, collaborative use of relays to 
form a beamforming system with the aid of perfect channel 
state information (CSI) and to provide physical-layer security 
is investigated. In particular, a decode-and-forward-based relay 
beamforming design subject to total and individual relay power 
constraints is studied with the goal of maximizing the secrecy 
rate. It is shown that the total power constraint leads to a 
closed-form solution. The design under individual relay power 
constraints is formulated as an optimization problem which is 
shown to be easily solved using two different approaches, namely 
semidefinite programming and second-order cone programming. 
Furthermore, a simplified and suboptimal technique which 
reduces the computation complexity under individual power 
constraints is presented. 

I. INTRODUCTION 

The open nature of wireless communications allows for the 
signals to be received by all users within the communication 
range. Thus, wireless communication is vulnerable to eaves- 
dropping. This problem was first studied from an information- 
theoretic perspective in [1] where Wyner proposed a wiretap 
channel model. Wyner showed that secure communication is 
possible without sharing a secret key if the eavesdropper's 
channel is a degraded version of the main chain. Later, 
Wyner's result was extended to the Gaussian channel in 
[3] and recently to fading channels in [4]. In addition to 
the single antenna case, secrecy in multi-antenna models is 
addressed in [5] and [6]. Liu et al. [7] presented inner and 
outer bounds on secrecy capacity regions for broadcast and 
interference channels. The secrecy capacity of multi-antenna 
broadcasting channel is obtained in [8]. Additionally, it's well 
known that that users can cooperate to form a distributed 
multi-antenna system by relaying. Cooperative relaying under 
secrecy constraints was recently studied in [10]— [12] . 

In [10], a decode-and-forward (DF) based cooperative pro- 
tocol is considered, and a beamforming system is designed 
for secrecy capacity maximization or transmit power min- 
imization. However, in this work, only total relay power 
constraint is imposed. In this paper, we extend the analysis to 
both total and individual constraints. Under individual power 
constraints, although analytical solutions are not available, we 
provide an optimization framework. We use the semidefinite 
relaxation (SDR) approach to approximate the problem as a 
convex semidefinite programming (SDP) problem which can 
be solved efficiently. We also provide an alternative method 
by formatting the original optimization problem as a convex 
second-order cone programming (SOCP) problem that can 
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Fig. 1. Channel Model 



be efficiently solved by interior point methods. Finally, we 
describe a simplified suboptimal beamformer design under 
individual power constraints. Simulation results are given to 
compare the performances of different beamforming methods. 

II. Channel Model 

We consider a communication channel with a source S, a 
destination D, an eavesdropper E, and M relays {R m }m=i 
as depicted in Figured] We assume that there is no direct 
link between S and D. We also assume that relays work 
synchronously by multiplying the signals to be transmitted 
with complex weights {w m } to produce a virtual beam point 
to destination and eavesdropper. We denote the channel fading 
coefficient between S and R m as g m <E C, the fading 
coefficient between R m and D as h m 6 C, and the fading 
coefficient between R m and E as z m E C. In this model, the 
source S tries to transmit confidential messages to D with the 
help of the relays while keeping the eavesdropper E ignorant 
of the information. 

It's obvious that our channel is a two-hop relay network. In 
the first hop, the source S transmits x s to relays with power 
-^[l^sl 2 ] = Ps- The received signal at R m is given by 



(1) 



where r\ m is the background noise that has a Gaussian distri- 
bution with zero mean and variance of N m . 

In the second hop, we employ decode and forward transmis- 
sion scheme. In this scheme, each relay R m first decodes the 
message x s and normalizes it as x' s — x s /P s . Subsequently, 
the normalized message is multiplied by the weight factor w m 



to generate the transmitted signal 



The output 



power of each relay R m is given by 



Mkl 2 l=* m x'J 2 l = \w„ 
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The received signals at the destination D and eavesdropper 
E are the superpositions of the signals transmitted from the 
relays. These signals can be expressed, respectively, as 
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where uq and n\ are the Gaussian background noise compo- 
nents at D and E, respectively, with zero mean and variance 
N . Additionally, we have defined h = [h* l: ....h* M ] T , z 



r * it 



] T , and w = [wi, ...wm] T where superscript 



denotes conjugate operation, and ( ) T and (-y denote the 
transpose and conjugate transpose, respectively, of a matrix 
or vector. The metrics of interest are the received SNR levels 
at D and E, which are given by 
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It has been proved that the secrecy rate R s over the channel 
between the relays and destination is 



R s = I(x s ;y d ) - I{x s ;y e ) 
= log(l+r d )-log(l + r e 
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where /(•; •) denotes the mutual information. In this paper, 
we address the joint optimization of {w m } with the aid 
perfect CSI and hence identify the optimum collaborative relay 
beamforming (CRB) direction that maximizes the secrecy rate 
given in ([9j. 

III. Optimal Beamforming under total power 

CONSTRAINTS 

In this section, we consider a total relay power constraint 
in the following form: ||w|| 2 = w^w < Pt- The optimization 
problem can now be formulated as follows: 



C s (h,z, Pt) — max log 
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where A max (A, B) is the largest generalized eigenvalue of the 
matrix pair (A, B). Hence, the maximum secrecy rate in (TT3b 
is achieved by the optimum beamforming weights given by 
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(14) 



where u is the eigenvector that corresponds to X max (NoI + 
Pyhh^j NqI + Ptzz^) and ? is chosen to ensure w„ pt w 0! , ( = 
Pt- 

Note that in the first-hop of the channel model, the maxi- 
mum rate we can achieve is 



(15) 



Ci - min log 1 + i£iii]_ 
m=l,...,M \ N„ 

Thus, the overall secrecy rate is 

Cdof, s = min(Ci,C s ) 
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Next, we provide some remarks on the performance of collab- 
orative relay beamforming in the high and low SNR regimes. 
For simplicity, we assume in the following that the noise 
variances at the destination and eavesdropper are iVo = 1. 



A. High SNR Regime 

In the high SNR scenario, where both P s , Pt — > oo and 

p, 
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lim = 1, we can easily see that 
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/,From the result in [6], we can see that 

lim C s = log A max (hh 1, , zz f ) 

P-r — >oo 



(18) 



= logP T + log(max|h t i/;| 2 ) (19) 

ip unit vector on the null space of z^. 

Based on these observation, at high SNR, we can choose 
the beamforming vectors to lie in the null spaces of the 
eavesdropper's channel vector. That is \^2 m —iZ m w m \ 2 = 
z'w = 0. In this case, the eavesdropper can not receive 
any data transmission from relays, the secrecy constraint is 
automatically guarantied. No secrecy coding are needed at 
relays. Now as Pt — > oo the weights optimization problem 
becomes 
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s.t z^w = 



log(P T ) + log(max |h^| s 
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since d2Tb is identical as ( fT9l , the null space beamforming 
is optimal at high SNR. Furthermore, the optimal null space 
beamforming vector can be obtained explicitly. 

Due to the null space constraint, we can write w = H^v, 
where denote the projection matrix onto the null space 
zl. Specifically, the column of Hjr are orthonormal vector 
which composed the basis of null spaces of z'. In our case, 
is a M * (M — 1) matrix. The power constraint w^w = 
vtH^H^v = vtv < P 
problem can be written as, 



Then the weights optimization 
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the optimum null space beamforming weights w is 
<;i is constant satisfied power constraint. 
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5. Low SNR Regime 

In the low SNR scenario , where both P s , Pt — > and 
lim = l.We can see that 
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Thus, in the low SNR regime, the direction of optimal CRB 
vector approaches the regular eigenvector to the largest regular 
eigenvalue of hh^ — zzL 

IV. Optimal Beamforming under Individual power 

CONSTRAINTS 

In a multiuser network such as the relay system we study in 
this paper, it becomes necessary to consider individual power 
constraint at the relays, which can be reflected by having 
| w m | 2 < p m Vm or |w| < p where we have | - | 2 to denote 
element-wise norm-square operation and p is a column vector 
containing the element {p m }. In what follows, the problem of 
interest will be maximization the term inside log function of 
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A. Semidefinite Relaxation (SDR) Approach 

We first consider a semidefinite program method similar to 
that in [13]. Using the definition X = ww', the optimization 
problem ( |3~T1 > can be write as 



N + Tr(hhtX) 

max ; ; — - 

x N + Tr(zrfX) 
s.t diag(3C) < p 

rank X = 1, and X >- 
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s.t ir(X(hh f - fzzt)) > N (t - 
diag(X) < p, 

rank X = 1, and X >- 



1). 



Where tr(-) represent the trace of a matrix and X >z 
means that X is a symmetric positive semi-definite matrix. 
The optimization problem in (l33l is not convex and may not 
be easily solved. Let us then ignore the rank constraint in 
(l33l . That is, using a semi-definite relaxation (SDR), we aim 
to solve the following optimization problem: 



max 
x.t 



t 



(34) 



s.t tr(X(hh t - izzt)) > N (t - 1), 
and diag(X.) < p, and X y 

If the matrix X opt obained by solving the optimization prob- 
lem in d34"i l happens to be rank one, the its principal component 
will be the optimal solution to the original problem. Note 
that the optimization problem in d34b is quasiconvex. In fact, 
for any value of t, the feasible set in ( f34b is convex. Let 
imax be the maximum value of t obtained by solving the 
optimization problem d34i >. If, for any given t, the convex 
feasibility problem 



find X 

such that tr(X(hh f - izz f )) > N (t - 
and diag(X.) < p, and X y 



(35) 
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is feasible, then we have t max > t. Conversely, if the convex 
feasibility optimization problem d35l l is not feasible, then we 
conclude t max < t. Therefore, we can check whether the 
optimal value t max of the quasiconvex optimization problem 
in (f34t is smaller than or greater than a given value t by 
solving the convex feasibility problem ( 1351 ). If the convex 
feasibility problem d3~5T l is feasible then we know t max > t. 
If the convex feasibility problem (f35T > is infeasible, then we 
know that i roax < t. Based on this observation, we can 
use a simple bisection algorithm to solve the quasiconvex 
optimization problem (f34b by solving a convex feasibility 
problem d35l)at each step. We assume that the problem is 
feasible, and start with an interval [I, u] known to contain 
the optimal value £ max . We then solve the convex feasibility 
problem at its midpoint t = (I + u)/2, to determine whether 
the optimal value is larger or smaller than t. We update the 
interval accordingly to obtain a new interval. That is, if t is 
feasible, then we set I = t, otherwise, we choose u — t and 
solve the convex feasibility problem again. This procedure is 
repeated until the width of the interval is smaller than the 
given threshold. Once the maximum feasible value for t max is 
obtained, one can solve 



minX tr(X) 
s.t ir(X(hh t - 
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and diag(~K) < p, and X >z 
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to get the solution X opt .([36]i is a convex problem which can 
be solved efficiently using interior-point based method. 

To solve the convex feasibility problem, one can use the 
well-studied interior point based methods. We use the well de- 
veloped interior point method based package SeDuMi [14]that 
produces a feasibility certificate if the problem is feasible and 
its popular interface Yalmip [15]. In semi-definite relaxation, 
the solution may not be rank one in general. Interestingly, in 
our extensive simulation results, we have never encountered a 
case where the solution X opt to the SDP problem has a rank 
higher than one. In fact, there is always a rank one optimal 
solution exists for our problem which we will explain later. 
Then we can obtain our optimal beamforming vector from 
principal component of the optimal solution X opt . 

B. Second-order Cone Program (SOCP) Approach 

The reason that this SDR method is optimal for above 
problem is that we can reformulate the problem as the second 
order cone problem. Thus, we have another way of solving it. 
The optimization problem(l30t is equivalent to 
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Observe that an arbitrary phase rotation can be added to 
the beamforming vector without affecting the constraint (38) . 
Thus, can be chosen to be real without loss of generality. 
Note (|38l can be write as 
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because h^w can be assume to be real, we may take the square 
root of the above equation. The constraint becomes a second- 
order cone constraint, which is convex. The optimization 
problem now becomes 



max t 



s.l \/-h t w> 



(40) 
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and 



|w| 2 <p 



As described in the SDR approach, the optimal solution of 
( |4TT i can be obtained by repeatedly check feasibility using a 
bisection search over t with the aid of well studied interior 
point method for second order cone program. Again, we 
use SeduMi together with Yalmip in our simulation. Once, 
the maximum feasible value i max is obtained, we can then 
solve following second order cone problem (SOCP) to obtain 
optimum beamforming vector. 



s.t 



-h f w > 



(41) 



and |w| < p 



Thus, we can get the secrecy capacity C s ,ind for the second 
hop relay beamforming system under individual power con- 
straint by the above two numerical optimization method. Then 
combined with first hop source relay link secrecy capacity 
C\. Secrecy capacity of decode and forward CRB system is 



Cdo/,1 



rnm(Ci, C s 



C. Simplified Suboptimal Design 

As was shown above, design the beamformer with indi- 
vidual relay power constraints requires an iterative procedure 
where, at each step, a convex feasibility problem is solved. 
We now try to get a suboptimal beamforming vector without 
significant computation complexity. 

we choose simplified beamformer asv = dw opt where w opt 



is given by ( f2Tb with 
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I 2 = Pt = I> 
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and we choose 
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where w op t,k,Pk is the fcth entry of w op t,p respectively and 

\w op tA 2 



k = are' max 

Ki<M 



Pi 



(43) 



Then,substitute this beamformer v into OTb to get the achiev- 
able suboptimal rate under individual power constraint. 

V. Simulation Results 

In our simulations we focus on the performance of 2nd 
hop secrecy capacity since the main concern of this paper 
is the design of the CRB system in second hop. We assume 
rayleigh distributed {g m }< {h m },{z m } is complex circularly 
symmetric gaussian with zero mean and variance ct 2 , errand 
a 2 z respectively. 

In Figj2] Fig[3] we plotted the second hop secrecy rate, 
which is the maximum secrecy rate our collaborative relay 
beamforming system can supported for both total and in- 
dividual relay power constraint. For individual relay power 
constraint, we consider the relays to have the same power 
budgets: Pi = Specifically, in FigfJJ we have ah = 3, 
a z = 1, Nq = 1 and Al = 5. In this case, the legitimate user 
has a stronger channel. In Figf3] we only change er/, = 1, 
a z = 2 that is the eavesdropper has a stronger channel. 
Our CRB system can achieve secrecy transmission even the 
eavesdropper's has a favorable channel condition. As can 
be seen from figures, the secrecy rate under total transmit 
power and that under individual relay power are very close 
to each other. And our two different optimization approach 
gives nearly the same result. It also can be seen that under 
individual power constraint, the simple suboptimal method 
suffers a constant loss as compared to SDR or SOCP based 
optimal value. 

In Fig|U we fix the relay total transmitting power as 
Pt = 10, and vary the number of collaborative relays. Other 
parameter is the same as Fig[3] We can see increasing the 
M will increasing the secrecy capacity for both total and 
individual power constraint. We also see that in some situation, 
increasing M will degrade the performance for our simplified 
suboptimal beamformer. 
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Fig. 2. 2nd hop secrecy rate vs. the total relay transmit power Pj< for 
different cases, eavesdropper has a weaker channel 
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Fig. 3. 2nd hop secrecy rate vs. the total relay transmit power Pt for 
different cases, eavesdropper has a stronger channel 
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Fig. 4. 2nd hop secrecy rate vs. Number of relays for different cases 
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